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The covariant spectator formalism is used to model the nucleon and the A(1232) as a system of 
three constituent quarks with their own electromagnetic structure. The definition of the "fixed-axis" 
polarization states for the diquark emitted from the initial state vertex and absorbed into the final 
state vertex is discussed. The helicity sum over those states is evaluated and seen to be covariant. 
Using this approach, all four electromagnetic form factors of the nucleon, together with the magnetic 
form factor, G* M , for the yN —» A transition, can be described using manifestly covariant nucleon 
and A wave functions with zero orbital angular momentum L, but a successful description of G* M 
near Q 2 = requires the addition of a pion cloud term not included in the class of valence quark 
models considered here. We also show that the pure 5-wave model gives electric, G* E , and coulomb, 
G* c , transition form factors that are identically zero, showing that these form factors are sensitive 
to wave function components with L > 0. 



I. INTRODUCTION 

Experimentally the main source of information on the 
internal structure of baryons lies in the electro- and 
photo-excitation of the nucleon, and is parametrized in 
terms of electromagnetic form factors. The precise elas- 
tic electron-protonpolarization transfer experiments un- 
dertaken at Jlab [2 H IU disclosed results at odds with 
previous data, triggering an intense discussion about the 
shape of the nucleon Jj, || (a review of the subject can 
be found in Refs. [H, of). Recent measurements of the 
7./V — > A transition raise further questions. How much 
of this process is due to the three valence quarks in the 
A, and how much to the meson-nucleon, or "pion cloud" 
contribution? Experimental progress has been impres- 
sive, with the recent accumulation of high precision data 
over a wide momentum transfer (Q 2 ) range. Just ten 
years ago, for example, the sign of the electric and mag- 
netic ratio for the jN — > A transition, was not even 
known beyond the photon point (Q 2 = 0). New pre- 
cise data have been collected from MAMI [8], LEGS Q, 

MIT-Bates [T3| and Jlab [UG3 in a re g ion Q 2 < 6 GeV2 
(q 2 = —Q 2 is the squared transferred momentum). 

This is the second in a series of papers using the covari- 
ant spectator theory to study the implications of mod- 
eling the baryon resonances as covariant bound states 
of three valence constituent quarks (CQ). In this model 
the structure of the CQ (including an anomalous mag- 
netic moment) arises from the dressed 7 — > qq interac- 
tions which give rise to the familiar meson structure of 
the vector dominance model, and these contributions are 
not included in the wave function of the nucleon. This 
language differs significantly from the light-cone formal- 
ism, where all of this vector dominance physics must be 
included in higher, non- valence components of the Fock- 



space expansion of the nucleon wave function (for further 
discussion, see [f|). In our first paper || (referred to as 
Ref. I) we showed that a very simple pure S-wave model 
of the nucleon, based on a covariant generalization of a 
simple non-relativistic SU(2) x SU(2) wave function for 
three valence CQ, could describe the four elastic nucleon 
form-factors very well. Our best model used 8 parame- 
ters: two to model the "radial" structure of the nucleon 
wave function, two to describe the anomalous magnetic 
moments of the up and down quarks, one to allow for an 
overall renormalization of the quark charge at very large 
Q 2 , and three to describe the vector dominance structure 
of the four quark form factors. Three of these parameters 
were fixed, leaving 5 to be varied during the fit. 

The present paper extends this model to the descrip- 
tion of the 7-/V — > A transition. As required by the model 
we use the same CQ form factors and nucleon wave func- 
tion as fixed in Ref. I. The only freedom remaining is 
the structure of the A wave function, and in this pa- 
per we study the consequences of the assumption that 
this wave function is a pure symmetric S-w&ve system 
of three valence CQ with total spin and isospin equal to 
3/2 (described with two parameters). We conclude that 
although the electromagnetic form factors of the nucleon 
alone may be described with such a simple ansatz based 
on orbital S'-waves only, the simultaneous description of 
all the 7iV — > A multipole transition form factors de- 
mands partial waves with L > 0. This is in agreement 
with the findings of the first non-relativistic quark mod- 
els which yield a magnetic dipole (Ml) form factor but 
no contributions to the electric (E2) and Coulomb (C2) 
quadrupole form factors [13] , since they did not include 
single quark D-states, either in the nucleon or in the A 
wave function. 

In Section II we introduce the A wave function, after 
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FIG. 1: Baryon-quark wave function defined in Eq. 12.lt . 

a short redefinition of the nucleoli wave function. In Sec- 
tion III the issues related to the basis for the polarization 
states are briefly reviewed. In Section IV the current 
corresponding to the N— > A electromagnetic transition 
is introduced, and in Section V the form factors are cal- 
culated. Section VI presents the results. Section VII 
summarizes the work and draws conclusions. 



II. RELATIVISTIC S-WAVE FUNCTIONS FOR 
BARYONS 

In thespectator framework 0, EI EE EE E3, El El, 
HE Hli [U , a general baryon with four- momentum 
P and mass M is described by a wave function for an 
off-shell quark and an on-mass-shell diquark-like cluster, 
defined by 

*(P,A;) = (m g -y 1 )- 1 (fe|r|P), (2.1) 

where T is the vertex function describing the coupling 
of an incoming on-shell baryon with mass M to an out- 
going off-shell quark and an on-shell quark pair (the di- 
quark). The quark pair is non- interacting with a contin- 
uous mass distribution varying from 2m q to infinity. All 
of the matrix elements in the spectator theory involve 
an integral over this mass distribution, and for simplic- 
ity this integral is approximated by fixing this mass at 
a mean value, m s (a parameter of the theory that scales 
out of the form factors, but can be fixed by deep in- 
elastic scattering, as described in Ref. I). The diquark 
four-momentum, k = P — pi, is constrained by its on- 
mass-shell condition k 2 = m 2 s . The quark has dressed 
mass rriq and four-momentum p\ . This wave function is 
shown diagrammatically in Fig. [T] 

Following Refs. I and [3], we write the baryon states 
in terms of the quark spin and "fixed-axis" diquark po- 
larization states, labeled by Ep. These vectors were dis- 
cussed in detail in Ref. I and [141 ] . We use them, instead 
of the diquark helicity vectors, since the latter depend 
not only on the magnitude of the diquark momentum, 
but also on its direction, while here we want to consider 
5-state orbital effects alone. Considering the "fixed-axis" 
polarization states we have shown [f| E3| that the wave 
function for the nucleon transforms as a Dirac spinor un- 
der a Lorentz operation. This implies that the model is 
covariant and assures the covariance of the electromag- 
netic current elements. The method was tested for the 



nucleon alone in Ref. I; here we test it for the A. 

A. S-wave nucleon wave function 

The manifestly covariant 5*-wave nucleon (with mass 
m) wave function was introduced in Ref. I 

$ NXn (P,k) = -±=^ N (P,k)$u(P,\ n ) 

- -^ N {P,k) ^ f P u{P,\ n ). (2.2) 

Here u(P, X n ) is a four-component Dirac spinor, and ipN 
is a scalar function that specifies the relative shape of 
both spin-isospin (0,0) and spin-isospin (1,1) diquark 
components. In the following we will sometimes sim- 
plify the notation by suppressing reference to the polar- 
ization A„ of the nucleon, and write the nucleon spinor 
as u(P, A n ) — > u(P). 

In Eq. (|2.2[) . (f>j gives the isospin states of the quark- 
diquark system (I = ±1/2 is the isospin projection of the 
nucleon) 

$ = e°v 

where £ l (i = 0, 1) represents the two diquark isospin 
states, and 

X + ^(J)=- X -i = (j)=p. (2.4) 

More details can be found in Ref. I. 

The spin-1 component of the wave function depends 
on the fixed-axis diquark polarization vectors e\ p = 
(eoi e x, £y, £z), with A = 0, ±1 the polarization index. 
The explicit expressions for e\ p are given in the next 
Section (Eqs. (|3 . 2[) for mjp = m). Here we emphasize 
that the polarization vectors are written in terms of the 
nucleon momentum P, instead of the diquark momentum 
k. This dependence gives <5 n the correct non relativistic 
limit and also assures that ty N satisfies the Dirac equa- 
tion 0. 

Note that (|2.2[) is written in terms of e* P , allowing the 
interpretation of $ n as an amplitude for an incoming 
nucleon and an outgoing diquark in the final state (see 
Fig. p. 

In this paper we only consider transitions from the nu- 
cleon to the A in which the diquark remains a spectator, 
with its spin and isospin unchanged during the transition. 
Since the diquark in the A must have spin-isospin quan- 
tum numbers (1,1) [as discussed below], only the spin- 
isospin (1,1) component of the nucleon wave function is 
needed. It is convenient to introduce a new notation for 
this component, and rewrite the (1,1) component of Eq. 
(gH) as 

= -4?ip N (P,k)(t>}sp*U a (p\ N ), (2.5) 



3 



where a is a vector index, and 



U a (P,X N ) = ^75 [7a - ~z: ) u(P,X N ) 



Ell 

m 



(2.6) 



Because ep ■ P = 0, this definition is equivalent to the 
one presented in [5[. However, the spinor (|2.6[) has the 
properties 



P a U a (P,X N ) = 
(to- P)U a (P,X N ) = 



(2.7) 



which are very convenient for the actual calculation of 
the matrix elements of the electromagnetic current. Note 
also that 



(m- P)$ NXN (P,k) =0. 



(2.8) 



The part of the wave function that depends on the 
magnitude of the relative momentum is modeled by the 
two parameter function used in Rcf. I: 



<MP,fc) 



A 



m s (Pi +xn){P2 + Xn) 



(2.9) 



Here No is a normalization constant and Pi with i = 
1,2 are range parameters in units of mm s . All of these 
parameters were fixed in Ref. I, and we use the same 
values in this paper. The dimcnsionless variable xn is 
defined as 



XN 



(to — m s ) 2 — (P — k) 2 



(2.10) 



Since the baryon and the diquark are both on-shell, the 
wave function (|2.9p can depend only on the variable 
(P — k) 2 (as required by the Hall- Wightman theorem) . In 
Appendix [G] we show that Eq. (|2.9[) assures the asymp- 
totic behavior for the nucleon form factors Ge and Gm 
will be 1/Q 4 times logarithmic corrections as expected 
from perturbative QCD (pQCD). 

This formalism is not restricted to the 5-wave case 
presented here. It can be extended to states with higher 
orbital angular momentum. This will be the subject of 
future work. 

The wave function (|2.5p has a very simple physical in- 
terpretation. It is a spin 1/2 state composed of spin 1 
and spin 1/2 constituents. This spin content is discussed 
in detail in Appendix [A] 



B. S-wave A wave function 

The iS-state wave function for the A (with mass M) 
is defined similarly to (1,1) component of the nucleon 
wave function in Eq. (|2.5p above. Non-relativistically a 
pure S'-wave spin 3/2 Delta state can be written as a 
direct product of spin 1/2 quark and a spin-1 diquark. 
Many details of the nonrelativistic construction of the A 



wave, and its relativistic generalization, are discussed in 
Appendix iBl In this section we summarize the results. 

In parallel to Eq. ()2.5j) . the S'-state Delta covariant 
wave-function can be written 



AA, 



XP,k) = -4>A(P,k)4> I ,e p *w f3 {P,X A ), (2.11) 



where 4>ii is the isospin part of the A wave function (in- 
cluding a diquark with isospin 1, and playing the same 
role as the nucleon isospin function </>}) with the isospin 
projections I' = ±1/2 or ±3/2, wp(P,X&) is the spin 
3/2 Rarita-Schwinger spinor- vector with spin projections 
Aa = ±1/2, ±3/2, and ?/>a(-P, &) is a scalar wave func- 
tion. In parallel with the nucleon definition (|2.5p . we 
define the wave function with a minus sign. For nota- 
tional simplicity, the spin indices of the diquark have been 
omitted from Eq. (|2.1ip implying that e" P — » e P with 
A = — 1,0, +1. The A wave function therefore consists 
of three components corresponding to the three different 
diquark polarizations. These polarizations are summed 
in the calculation of the transition form factors, as dis- 
cussed below. The diquark spin polarization vector is, as 
in the nucleon case, a function of the A mass and mo- 
mentum. It will be given explicitly in Eqs. (|3.2[) with 
m H = M. 

Since wp(P, Aa) are Rarita-Schwinger spinor- vectors, 
they satisfy the usual constraint conditions [2J, |25| 



(M— P)wp(P,Xa) = 
P^w (P,X A ) = 
1 P w p {P,Xa) = 0. 



(2.12) 



Therefore, the wave function (|2.1ip satisfies the Dirac 
equation 



(M- P)V A x & (P,k) = 0, 



(2.13) 



showing that ^aaa (Pi k) has no lower (negative energy) 
components in its rest system. 

The isospin wave function (f>j, can be written as 



h> = (t ■ e*) x 1 ' 



(2.14) 



where £ is the isospin vector of the diquark (identical 
to the one used for the nucleon), x 1 is a 4 x 1 isospin 
state, and T % denotes the 2x4 matrix corresponding to 
the 3/2 — > 1/2 isospin transition operator. The specific 
form of these operators is given in Appendix [B] 

As in the nucleon case, tpA can be expressed as a func- 
tion of 



XA 



(M - m s f - (P - k) 2 
Mm* 



In particular, we use phenomenological ansatz 

Ni 



^ A (P,k) 



m s (a 1 + xa)(«2 + Xa) 2 ' 



(2.15) 



(2.16) 



where on (i = 1,2) are range parameters in units of 
Mm s and N± a normalization constant. We note that 
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the power of xa in the denominator differs from the cor- 
responding one in the phenomenological ansatz for the 
nucleon case. In Appendix [U] we show that this choice 
for the A wave function, together with the nucleon wave 
function (|2.9[) . give the expected 1/Q 4 pQCD limit for 
the dominant form factors of the electromagnetic N — > A 
transition. 



C. Covariant spin projection operators 

Both the nucleon and A wave functions have a generic 
structure 

*hx h = <f> H (P,k) e P *V Ha {P,\ H ) (2.17) 
where H = N, A, and 




U a (P,XN) nucleon 
w a {P,\ A ) A 

h=t(j N (P,k)(j) 1 I nucleon 



V2 



-l/>A{P,k)(f>i 



A. 



(2.18) 



Furthermore, in both cases 



P a V a (P,X H ) = 0- 



(2-19) 



These similarities allow us to make some interesting gen- 
eral statements about the nucleon and A wave functions. 

The condition (|2.19[) means that, in the rest frame of 
the hadron, the vector V a has spatial components only. 
In this subspace the identity operator is 



0=9 0=9 13 



73 



(2.20) 



where we use the notation = m and wia = M. This 
subspace is spanned by two projection operators: 



y a0 
1 1/2 

where, using the notation 



'3/2 y ■ 



ppo 



m 



(2.21) 



(2.22) 



H 



the two operators are 



V 



a/3 
1/2 
-, a 
3/2 



n%{p) 

P$(P) = u 



a/3 



1 ^0. 



(2.23) 



In Appcndix[C]we show that these operators are relativis- 
tic generalizations of the spin 1/2 and spin 3/2 projec- 
tion operators (for a particle of mass rn#). They op- 
erate only in the 3x3 subspace of space-like vectors 
~a _ y a _ pa (p. v y m ^ _ These operators are well known 
in the literature, in others contexts (the operator V1/2 is 
sometimes denoted V\\) [26l 



As expected, the nucleon and A wave functions are 
eigenvectors of the spin 1/2 and spin 3/2 operators: 



K/2 w p 



u a 
= 



P$U = O 

^3/2 W(3=W 



(2.24) 



and the orthogonality of these projection operators 



i aP v. 



0U 

3/2 







(2.25) 



implies that the two wave functions are orthogonal, as 
expected. We will see later that a generalization of this 
condition is useful in proving current conservation in the 
electromagnetic N — ► A transition process. 

In work already underway [28| these operators will also 
play an important role in the extension of the formalism 
to D-wave states, for both the nucleon and the A. 



III. FIXED-AXIS DIQUARK POLARIZATIONS 

A. Definition of the state vectors 

The spin-1 fixed- axis polarization vectors were intro- 
duced in Ref. I, and some of their properties discussed 
and derived in [14j |. (These are really axial vectors, 
but for simplicity we will drop the word "axial" in the 
subsequent discussion.) These vectors are denoted s^ p , 
where A = 0, ±1 is the spin projection in the direction of 
the baryon three-momentum, P, with P — {£ p ,P} the 
baryon four-momentum and £ p — \J m 2 H + P 2 the baryon 
energy. In the baryon rest frame the fixed axis may be 
chosen to be in any direction. Choosing the z direction 
the polarization vectors are 



■OPo 
±Po 



(0,0,0,1) 
^(0,TF-*,0). 



(3.1) 



where the baryon four-momentum in its rest frame is 
denoted Pq = {m#, 0,0,0}. These polarization vectors, 
when used in Eqs. (|2.5|) and (|2.11j) . give the correct non- 
relativistic limit for the nucleon and A wave function. 

To write the baryon wave function in a frame where 
the baryon is moving, a boost in the direction of motion 
(^-direction by convention) of the baryon is needed. For 
this choice, P = (5p, 0, 0, P), and the polarization vectors 
become 



OP 



±p 



— (P,0,0,£ P ) 
7f(0,T,-*,0) , 



and satisfy 



£* XP ■ E\'p — —8\\> 



sxp-P = 0. 



(3.2) 



(3.3) 



In the following we will use the variable P- for the 
momentum in the initial state and P + for the momentum 
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in the final state. Also, the initial state will be a nucleon 
and the final state will refer either to nucleon or to a A. 

In the Breit frame, for a transition from an initial state 
of mass m to a final state of mass M, with a momentum 
transfer q = P + — P_ , we write 



£+,0,0,-q L ) , 



(3.4) 



where E + = J M 2 + \q\ and E. 



ql = Q 2 



(M 2 - m 2 ) 



= \ m 



2 \2 



\q\ with 



2(M 2 + m 2 ) + t 



(3.5) 



Eqs. (|3.4p hold for both equal masses M = m (q^ = 
\J Q 2 ) and unequal masses. 

In the Breit frame, according to Eqs. (|3.2p . we have for 
the initial state 



1 



-0P_ 

-±p_ 



V2 



(0,Tl,-i,0), 



(3.6) 



and for the final state 

fc 0P+ " 



M(hL,0,0,E, 



-±p 



V2 



(3.7) 



Note that the polarization vectors e P± from Eqs. (|3.6[) - 
(|3.7[) refer to the Breit frame only. 

Starting from the Breit frame we can then change to 
an arbitrary frame by means of a suitable Lorentz trans- 
formation A. The details of this transformation are dis- 



cussed in Ref. 14| . The Breit frame momentum P^ (with 
P = P±) is then transformed into P ,fi according to 



(3.8) 



Due to the four- vector character of e P , the polarization 
vectors in the new frame Em, parametrized by A, become 



'XPi 



A"„e 



(3.9) 



for each polarization A. In this notation only the momen- 
tum index distinguishes the arbitrary frame (momentum 
P') from the Breit frame (momentum P) polarization 
vector. 

The same transformation A acts on both final and ini- 
tial momenta P±. The transformation (|3.9p combined 
with the transformation law of the Dirac spinors and the 
Rarita-Schwinger states (see Appendix B) implies that 
the baryon wave functions transform as Dirac spinors. 
The demonstration of the same property for the A wave 
function follows the lines of the presented in the Ref. [l4[ • 

Fixed-axis polarization vectors Ep are different from 
the helicity vectors used in Ref. 29], which we denote 



here by 77. The latter depend on the diquark momentum 
k, and therefore on its direction, satisfying 77 • k — 0. The 
helicity vectors 77 can be related to our fixed-axis polariza- 
tion states by a rotation In fact, with an appropriate 
redefinition of the vertex function T, a wave function us- 
ing fixed-axis states can be made exactly equivalent to 
another wave function using helicity states. In the case 
of an initially totally spherical symmetric wave function, 
the transformation from diquark fixed-axis polarization 
states to direction-dependent or helicity states gives a 
vertex function accompanying the helicity vectors just 
the right angular dependence on the diquark momentum 
to cancel the dependence introduced by the helicity states 
77 [14| . Conversely, a spherically symmetric vertex func- 
tion r, like the one used here, if taken together with 
the direction-dependent diquark helicity states 77, would 
result in a wave function without spherical symmetry. 
Since here we want to investigate the consequences of 
spherical symmetric wave functions only, it becomes nat- 
ural to write these wave functions in terms of fixed-axis 
diquark polarizations. 



B. Importance of the collinearity condition 

We emphasize that matrix elements of states that in- 
clude fixed-axis polarization vectors must first be con- 
structed in a frame in which the incoming and outgoing 
states have collinear three-momenta, and only after this 
has been done can the matrix elements be transformed to 
an arbitrary frame. If matrix elements are constructed 
in this order, they will be both unique and covariant, 
but if they are not constructed in this order, they will 
be neither covariant nor unique. A simple example of 
the problems encountered if one does not start with a 
collinear frame is developed in Appendix [Dj Our failure 
to emphasize this point in our original presentation of 
these ideas lead to a criticism of Kvinikhidze and Miller 
[30j ] , which we addressed completely in Ref. [l4| • 

Why is a collinear frame required? We will see in the 
next section that the matrix elements we calculate as- 
sume that the diquark is a spectator which does not par- 
ticipate in the interaction. Hence, that the polarization 
of the diquark emitted by the initial baryon must be the 
same as the polarization of the diquark absorbed by the 
final baryon. There are not two distinct diquarks, but 
one and only one, diquark. Therefore, using fixed-axis 
polarization states, we must be certain that the polariza- 
tion states of the diquark emitted from the initial ver- 
tex, and of the diquark absorbed into the final vertex, 
are defined with respect to the same axis, and only in 
the collinear frame are we certain that the definitions 
of the polarization of the incoming and outgoing diquark 
(the same diquark) are consistent with a single direction. 
Therefore, if we happen to be presented with an interac- 
tion in which the initial baryon three-momentum P_ is 
not parallel to the final baryon momentum, P + , we must 
first transform the matrix clement of the collinear frame, 
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construct the matrix element, and then transform back 
to the original frame. 

Fortunately, given a any initial and final momentum 
configuration, there is always a Lorentz transformation 
that will boost and rotate both states to a collinear frame 
with the three-momenta P± in the z direction, so the 
need to define the states in collinear frame imposes no 
limitation. Our construction is similar to the definition 
of two-particle helicity states in the two-body center-of- 
mass by Jacob and Wick (3lj . 



P4 



N, A. 



/ 



k 



N 



FIG. 2: Relativistic impulse approximation. 



IV. MATRIX ELEMENTS OF THE CURRENT 



A. Definition of the current 



where r 3 is the quark isospin operator. In this work we 
adopted the quark form factors from Ref. [B[ 



Consider the electromagnetic transition from an initial 
state *S?i (mass m) and a final state (mass M). The 
relativistic impulse approximation (RIA) to the transi- 
tion current in the spectator formalism, shown in Fig. [21 



J%(P+,P-) = */(P+,fe)i^i(P-,fc). (4.1) 

This is covariant, but we will study it in the collinear 
Breit frame where the fixed-axis polarizations are con- 
sistently defined, as discussed above, and the incoming 
and outgoing momentum, P T , are already defined in 
Eq. (33J. 

In the RIA, the photon couples to each quark through 
the diagram shown in Fig. O The factor of 3 in (|4.1|) 
comes from isospin invariance, which allows us to express 
the sum of the three diagrams in terms of a single inte- 
gral multiplied by 3. All intermediate states are taken 
into account by summing over the diquark spin-1 polar- 
izations and integrating over all positive on-mass-shell 
diquark (spectator) states with energy E s , using 

^ (4-2) 



There is, in principal, an integration over all spectator 
masses, m s , but this integral is replaced by the value of 
the integrand at some (unknown) mean value m s , which 
becomes a parameter of the model. 

The quark current jj (which is isospin dependent) is 
decomposed into its Dirac and Pauli terms, 



% = hi* 



■32- 



2m 



(4.3) 



where ji [i = 1, 2) are the quark form factors, defined in 
Ref. I. 

These form factors j\ and j% include the quark struc- 
ture and can be decomposed into isoscalar and isovector 
components: 

ji = \h+{Q 2 ) + \h-{Q 2 )T3 

32 = \h+{Q 2 ) + \f2-{Q 2 )r^ (4.4) 



h±(Q 2 ) 
f2±(Q 2 ) 



= A + 



(1-A) 
1 + Q 2 /m 2 v 

d± 

1 + Q 2 /ml 



c±Q 2 /Mj 
(l + Q 2 /M 2 ) 
(1 ~ d±) 



1 



(4.5) 



In these expressions m v and Mh are vector meson masses 
that represent the dominant contributions from the vec- 
tor dominance model (VDM). The lower mass, m v = m p 
(or m w ), describes of the two pion resonance (three pion 
resonance) effect and Mh , fixed as 2m (twice the nucleon 
mass), takes account of all the large mass resonances. 

The parameter A is fixed by the deep inelastic scat- 
tering (DIS) distribution amplitudes and can interpreted 
physically as a scaling of the quark charges in DIS limit. 
All the other parameters are presented in Table HI They 
correspond to a previous application of the covariant 
Spectator theory to the description of the nucleon form 
factor data with only an 5-state in the nucleon wave func- 
tion [5j . We know from the beginning that the restriction 
to orbital S- waves is a considerable simplification, baryon 
ground state to its first resonance, but it is interesting to 
see exactly what are the consequences of such a simple 
assumption. 

Note that all of these parameters are fixed by the nu- 
cleon data: elastic form factors and DIS. The quark elec- 
tric form factors f\± are normalized to 1 at Q 2 = 0, in 
order to reproduce the quark and nucleon charge. The 
quark magnetic form factors are normalized by proton 
and neutron magnetic moments /2+(0) = k + = 1.639 
and /2-(0) = «- = 1.823; see Ref. 5] for a detailed 
discussion. 

Due to the relation T\ti — only the isovector com- 
ponents of the current contributes to the — > A tran- 
sition. Then, in the following discussion we need the 
isovector current only, which can be written 



3i 



T3 



i-r + h- 



2m 



73 
2 ' 



(4.6) 



with the isovector quark form factors /i_ and fa- defined 
as above. 
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Model /3i, /?2 c+,c_ d+,d- X,m s /m Nq,x 2 



I 


0.057 


2.06 


-0.444 


1.22 


10.87 




0.654 


2.06 


-0.444 


0.88 


9.26 


II 


0.049 


4.16 


-0.686 


0.547 


11.27 




0.717 


1.16 


-0.686 


0.87 


1.36 



TABLE I: Parameters of the nucleon wave function (/3i,/32 
and iVo) and quark form factors. In each case we kept 
k+ = 1.639 and k_ = 1.823 in order to reproduce the nu- 
cleon magnetic moments exactly. The difference in the \ 2 
is mainly due to the description of the neutron electric form 
factor. Model I preserves the isospin symmetry for the quark 
electric form factor /i+ = but cannot describe neutron 
electric form factor data. See details in Ref. [jj. 



B. Diquark polarization sum 

The next step in the general reduction of the transition 
current is to carry out the sum over the diquark polar- 
izations. Using the generic notation of Eq. (|2.17| . the 
current is written 



j£(p + ,p_)=; 



P (P+,fc)T 3 i (P-,A;) 



D f3 a 



xV f p(P + ,\ + )j!tV ia (P-,\-), (4.7) 



where we assume that only the isovector quark current 
contributes to form factor (true for the 7iV — > A tran- 
sition) and have written the diquark spin sum as the 
operator 



t \P + t XP- 



(4.8) 



that is evaluated in Ref. [lH and Appendix [El The final 
result shows that depends only on the momenta and 
masses of the two states, and can be written 



D'"' = - | g> 

a [ P^ 



p»pv 

' P+ P- 
P+ P 



(4.9) 



M 2 



-PI 



PI 



P+-P- 



P v _ 



where the factor a is 



a 



Mm 



P + ■ P_ [Mni + P + ■ P_ 
Note that D^ v satisfies the conditions 

P +t ,D^ = D^P- U = . 



(4.10) 



(4.11) 



current (from the quark anomalous magnetic moments). 
The Pauli current is always conserved, independent of 
the asymptotic states considered. To reduce the Dirac 
current we use the facts that the initial and final states 
both satisfy the Dirac equation, and that the charge form 
factors of the quark depend on q 2 and can be factored out 
of the integral 



\ J k 



= §(M-m)/i_ / VfT 3 ^ l (4.12) 

If the masses are equal, the condition is automatically 
satisfied, but for unequal, masses the states must be or- 
thogonal 



E / */ r 3^ 



0. 



(4.13) 



We can also write Eq. (|4.12|) using the notation of 
Eq. (14~T1) 



,{P + ,k)TZ<t>i{P-,k) 



x V ffi (P+,k + )D Pa V ia (P-,X-). (4.14) 

For the — > A transition, we can use the projection 
operators to prove orthogonality. Using the fact that the 
N and A states are eigenvectors of the spin- 1/2 and spin- 
3/2 projection operators, we can write 



V A p(P + ,\ + )DP a V Na (P^,\-) 



(4.15) 



= V A ^(P+, A+) |P$(P+) D Pa P«J 2 {P.)\ V Nv (P-,X-) , 

where V 3 / 2 (P+) and Vi/ 2 (P-) are the projection opera- 
tors of Eq. f2"^3"|) with P P+ and P -» P_, respec- 
tively. We show in Appendix [C] that the operator in 
square brackets is zero: 



Ptf 3 (P+)D Pa P1%(P-)=0 



3/2 



/2 



(4.16) 



This is the generalization of the orthogonality relation 
(|2.25[) and proves the orthogonality of the wave functions 
for all momentum transfers, q. Due to this orthogonal- 
ity between the initial and the final states, the additional 
-fiq^/q 2 term used in the definition of the current in Ref. I 
vanishes in this application. This is why we did not in- 
clude that extra term in Eq. 



V. 7 AT — > A TRANSITION 



Current conservation 



Current conservation requires that q^J^ — 0. To see 
if this condition is satisfied, we consider separately the 
Dirac current (from the quark charges) and the Pauli 



We will now apply the formalism of the previous sec- 
tions to the study of the electromagnetic iVA transition 
which has a simple interpretation in terms of valence 
quark structure: the A is a result of a spin flip of a sin- 
gle quark in the nucleon. It is then understandable that 
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the magnetic dipole multipole Ml dominates the transi- 
tion for low Q 2 , while the electric E2 and the Coulomb 
C2 quadrupoles give contributions of only a few percent. 
For large Q 2 however, according to perturbative QCD es- 
timations, we expect equal contributions from Ml and E2 
[32| . At the present the Q 2 scale of the pQCD regime is 
not known, which motivates calculations within models. 



A. Simplification of the transition current 

The transition current (|4. 7[) can now be simplified. Us- 
ing the notation of Eqs. (|2.18p and (|4.9[) and changing 
the integration variable from k to k/ m s we obtain an in- 
tegral independent of the diquark mass m s (due to the 
wave functions normalization which goes with l/m s ) and 
the diquark energy factor E s /m s = y/l + k 2 /m 2 (which 
altogether cancel the w? s dependence of the element d 3 k). 
Factoring out the isospin factors gives 



where 



P-) = - 5$5(0i')W} / i>A(P+,k)i> N (P-,k) 



xw (P+,X+)j!:U a (P.,X-)D 



(3a 



(5.1) 

where j£ is the defined in terms of the isovector part of 
the quark current (|4.6|) : 



# = f 



1-V + /2 



2m 



(5.2) 



which is the only part of the quark current to contribute 
to the transition amplitude. The isospin matrix element 
is evaluated using the properties of the isospin matrix 
transition T % (between spin 1/2 states and 3/2 states). 
Summing over the isospin projections mi of the diquark 
isospin vector gives 



2V2 , 



3//' 



(5.3) 



Next, using the fact that the initial and final states both 
satisfy the Dirac equation, we may reduce the Pauli form 
of the current using the Gordon decomposition 



M 



2m ^ \ 2m / m ^ ^ 

where P' 1 is the average of the initial and final momen- 
tum, defined in Eq. (|5.12p . We already saw in the discus- 
sion of gauge invariance [leading up to the identity (|4.16j) ] 
that the matrix element of the identity operator is zero, 
and hence the P^ term does not contribute. This allows 
us to collect the quark charge and anomalous magnetic 
moment contributions into a single term, giving finally 



j£ N (P + ,P_) = 6 ri f v I ip A (P +) k)i> N (P-,k) 



x wp (P+ , A+ ) r U a (P_ , A_ )D? a , 

(5.5) 



fv = h 



M 



2m 



-h 



(5.6) 



is a particular linear combination of the quark from fac- 
tors that can be factored out of the integral because it 
depends on Q 2 only. 

Equation (|5.5|) includes the explicit conservation of the 
z-projection of the isospin. This means that the model 
predicts that the amplitude is the same for both isospin 
channels: 7*p — *■ A + and 7*n — > A . 

The nucleon wave function ip^ is normalized to one, as 
required by the charge conservation at Q 2 = @. Sim- 
ilarly, also the A wave function (|2.11[) is constrained, in 
the rest frame where P = (M, 0, 0, 0), by the charge con- 
dition (excluding the isospin states from the wave func- 
tion): 

Ql = [ *A(P,*)il*A(P,*) 



1 + T 3 



\^A(P,k)\' 



(5.7) 



where j\ = \ + \t%. The isospin operator T3 is defined 



Ta = *Y, T i 7 * Ti = 

i 

Equation (|5.7p gives the correct A charge if 



3 

10 

0-10 

-3 



(5.8) 



|^A(P,fe)| S 



(5.9) 



This condition determines also the normalization con- 
stant for the wave function in Eq. (|2.16D . 



B. Form factors: Generalities 

The N — > A electromagnetic transition current (ex- 
cluding the electron charge e and ignoring the polariza- 
tions of the nucleon and the A) is given by 

J» = w p (P+)T^(P, qh 5 u(P-), (5.10) 

where the general form of the transition vertex is 

T^(P, q) = gV Gi + <fP»G 2 + <zVG 3 - g^G 4 . 

(5.11) 

The variables P and q are respectively the average of 
baryon momenta and the photon momentum: 



p = ±(p+ + P-) 

q = P+-P- 



(5.12) 
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The form factors Gi, i — 1, ..,4 depend only on Q 2 = 
—q 2 - Due to current conservation, q^T^^ — 0, only three 
of the four form factors are independent. In particular, 
we can write G4 in terms of the first three form factors 



Ga 



(M + m)Gi + i(M 2 -m 2 )G 2 



Q 2 G 3 , 



(5.13) 



and adopt the structure originally proposed by Jones and 
Scadron [33j . 

The parametrization (|5.1ip is not directly comparable 
to experimental data. More convenient for that purpose 
are the magnetic dipole (M) , electric quadrupole (E) and 
Coulomb quadrupole (C) form factors defined as 



G* M (Q 2 ) = k { [(3M + m)(M + m) + Q 2 



Gi 
M 



g* e (Q 2 



+ {M 2 -m 2 )G 2 ~2Q 2 G 3 } 
K {(M 2 -m 2 -Q 2 )^ 

+ (M 2 -m 2 )G 2 -2Q 2 G 3 } 
k |4MGi + (3M 2 + to 2 + Q 2 )G 2 

+ 2(M 2 -m 2 -Q 2 )G 3 }, 



where 



3(M + to) 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



The three form factors G* (a = M, E, C) are related to 
the magnetic, electric, and Coulomb (or scalar) multipole 
transitions, respectively. 



C. Form Factors: Application of the model 

Substituting for Up in (|5.5j) . and suppressing the 
isospin conservation factor Sjij, gives immediately 

J^(P +J P_) =wp(P + ,X+)0^ 5 u(P-,X-), (5.18) 

where 

= j-f v J [^ A (P + ,fc)V7 W (P_,fc)] 7 M C /3a 7a(5.19) 

This is easily reduced; the work is given in Appendix 
IFl The final results for the form factors of a transition 
between S'-wave nucleon and A states are 



C-(Q 2 ) = G^(Q 2 ) = 0, 



where 



(5.20) 
(5.21) 

(5.22) 



Note that I is the only factor which depends on the scalar 
wave functions and it is Lorentz scalar (frame indepen- 
dent). 



VI. RESULTS 

Before we present the numerical results we focus on 
the analytical structure of the results (|5.20|) - (|5.21|) . Note 
that the electric and Coulomb quadrupole form factors 
vanish in this calculation, which is restricted to S-wave 
orbital states. This result is consistent with quark mod- 
els based on quark S'-wave states [13j. According to the 
literature, the presence of multipoles E2 and C2 is a sig- 
nature of the nucleon and/or A deformation. Our result 
of identically vanishing quadrupole form factors is the 
consequence of considering only S-states for both nucleon 
and A wave functions, and consequently a nucleon and a 
A with spherical form. The inclusion of higher orbital 
momentum components would generate non-vanishing 
electric and Coulomb quadrupoles, as we will confirm 
in forthcoming work. It is worth noticing that the ex- 
perimental results for the transition multipoles indicate 
a contribution of E2 and C2 of the order of a few percent 
at Q 2 = 0, consistent with a small angular momentum 
component in the wave function. 

Next, look at the magnetic form factor, G* M , at Q 2 = 0. 
Substituting for f v gives 



Gm(0) 



2 m 



3^ \ M + m 



+ k_ 1(0). (6.1) 



The isovector magnetic moment K- was fixed by the nu- 
cleon magnetic moment in Ref. § 



K- = f(// p -//„)- 1 = 1.823, 



(6.2) 



For future discussion we write Eq. (|6.1[) as 



G* M (0) = 



V3 



fjjp fJ, n 



lM-m 



3 M + m 



1(0). (6.3) 



When the experimental nucleon magnetic moment is used 
in (16.11) one has 



where 



G* M (0) = 2.071(0). 



1(0) = / ^aVw 

Jk 



(6.4) 



(6.5) 



Note, however, that we are working with normalized wave 
functions 



= 1, 



Q 2 =0 



= 1. 



> 2 =o 



(6.6) 



Because of this conditions, the integral 1(0) is limited, in 
its absolute value, by the Holder inequality, the version 
of the Cauchy-Schwartz inequality for an integral 



V'aV'. 



N 




(6.7) 



10 



In particular for Q 2 — we obtain 

|X(0)| < 1. (6.8) 

Choosing the positive sign (consistent with our defini- 
tions of the scalar wave functions) we conclude that 

G* M (0) < 2.07. (6.9) 

Since the experimental value is considerably larger than 
this limit, 

G* M (0) =3.02 ±0.03, 

we see that the Spectator quark model, in impulse ap- 
proximation, can, at best, only describe 69% of the 
77V — > A transition form factor G* M at the photon point. 

This underestimation of G* M (0) is an universal prop- 
erty of all constituent quark models. It has been previ- 
ously reported in the literature [H [H, HE 113, S H, H| • 
Naive 5- wave non-relativistic quark models with SU (6) 
symmetry (and no dynamical effects included) [l3|, pre- 
dict 

2^2 



G* M (0) 



2.30. 



(6.10) 



Including kinematic effects, Ref. [4l| obtains 



G* M (0) 



2V2^2E N (m + E N ) 



M 



2.04. 



.11) 



where En is the nucleon energy at threshold in the A 
rest frame. For a review of the constituent quark models 
predictions see Ref. [l3j . 

Our model for a quark-diquark system differs from 
other quark models; our quarks are not static, as in Eq. 
(|6.10p . and our magnetic form factor is related to both 
Hp and not only to [L v as in Eq. (|6 . . Neverthe- 
less, comparing Eqs. (|6.3p and (|6.1ip we can see that 
very different descriptions can lead to the similar results 
if the same constraints are considered (normalization of 
the wave functions). For completeness we add that cal- 
culations based on QCD sum rules also lead to an under- 
estimation of G* M (0) although these models do not apply 
to the Q 2 = region [42[ ■ Similar results are obtained us- 
ing Generalized Parton Distributions (GDP) [H El IZ5I] . 
These models extrapolate the Parton Distribution Func- 
tions from Deep Inelastic Scattering to intermediate en- 
ergies. For low Q 2 G* M is underestimated by 20-30% 

mm. 

The failure of quark models to describe the 7^ — > A 
transition at threshold shows their limitations, which 
stem from taking constituent quarks as the only relevant 
degrees of freedom. Quark wave functions can be normal- 
ized to correctly describe nucleon and A static charges, 
but fail in the description of the dynamical 7iV — > A 
transition which does not involve a charge density, but 
instead a transition charge density. The magnitude of 
the difference between the quark model result, which we 
label the Bare result, and the experimental result, may 
be due to pion field contributions, and is a manifesta- 
tion of the strong correlation between the A and the 7rN 
system. 



O.N 



0.6^ 



0.4 



0.2 - 



Model 1(B) 

Model 11(B) 



"f- 

I 



012345678 
Q 2 (GeV 2 ) 

FIG. 3: Result of the fit the A wave function parameters 
to the data Q 2 > 2.9 GeV 2 (where G* M is expected to be 
very small). The nucleon and quark parameters are given by 
model I and II (see Table HJ. Delta parameters are presented 
in Table HU Data from CLAS/Jlab [UIH, DESY [H and 
SLAC 541. 



Decomposition into Bare and Pion Cloud form 
factors 



Following the previous discussion, we decompose G* M 
into two contributions: 



Gm(Q 2 HGm(Q 2 ) + gmQ 2 



(6.12) 



The term Gf f is the Bare form factor: the contribution 
of the quark core given by the quark model under con- 
sideration. The term G\,j is a contribution due to the 
pion field: the contribution from any diagram involving 
a photon and pion loops. Our spectator quark model can 
predict G^ only. 

There are two kinds of descriptions that take into ac- 
count the effects of the pion field explicitly: dynamical 
models and low momentum Effective Field Theories or 
Chiral Perturbation Theories. Here we focus on dynami- 
cal models 0, [49|, [5(| because these models can be used 
to describe the entire momentum region over which data 
is available. 

A dynamical model uses hadronic degrees of free- 
dom and a coupled channel method to derive transi- 
tion amplitudes involving initial and final meson-baryon 
and photon-baryon states. The transition amplitude 
can be decomposed in two components: (i) the back- 
ground or non-resonant amplitude which is the solu- 
tion of an Lippmann-Schwinger-like equation with a non- 
resonant interaction kernel; and (ii) the resonant ampli- 
tude which includes the contributions from dressed in- 
termediate baryon resonance states. The non-resonant 
interaction kernel that generates the non-resonant back- 
ground includes direct couplings of the photon or mesons 
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with the baryons, and may also include meson rescatter- 
ing described by intermediate vector mesons. The res- 
onant part is generated by dressing the s-channel pole 
terms generated by vertex functions describing the cou- 
plings of photons or mesons to the bare baryon pole. 
Both the non-resonant direct coupling terms and the res- 
onance vertex functions are parametrized by simple phe- 
nomenological expressions with parameters adjusted to 
fit the pion-nucleon and pion photo- production data. A 
review can be found in Refs. (l3l. l46l |47||. 

The effective contribution of the pion cloud depends on 
the particular model. The models of Sato and Lee (SL) 
[H , Dubna-Mainz- Taipei (DMT) [H predict that the 
pion cloud will give an important contribution at Q 2 = 
that falls quickly with increasing Q 2 . The Utrecht-Ohio 
model [50] in opposition predicts small contributions for 
Q 2 and more significant contributions of Q 2 ~ 2 
GeV 2 for the pion cloud. According to Ref. [51], pion 
cloud effects give 33% of the total contribution at Q 2 = 0, 
and less than 10% for Q 2 > 4 GeV 2 . 

Pure quark models include no pion cloud effects and 
can give a complete description only at higher Q 2 where 
contributions from the pion cloud are negligible. Sup- 
porting this interpretation, our numerical calculations 
show that we can fit the region Q 2 > 2.5 GeV 2 , but not 
the low Q 2 region. A fit to the higher Q 2 data (Q 2 > 2.9 
GeV 2 ) is presented in Fig. H The A wave function pa- 
rameters obtained from the two fits shown, referred to as 
Models 1(B) and 11(B), are given in Table[TIl We omitted 
from this Table the parameters of the nucleon wave func- 
tion model which also enters the calculation, since that 
wave function was already fixed by the nucleon form fac- 
tors and DIS results [f| , and those parameters were shown 
already in Table |TJ 

From the figure we conclude that with no explicit pion 
cloud we can explain about 55% of G* M at Q 2 = 0. This 
contribution is lower that the upper limit of Eq. (|6.4[) . 
For each model, the theoretical quantity T(0) is a mea- 
sure of the extent to which a model approaches its the- 
oretical upper limit, and the deviation of G^(0)/3 from 
the experimental value 1 is a measure of the quality of 
the fit to the data at Q 2 — 0. 

To compare our model with the data over the entire 
range of Q 2 , we need a parametrization for GJ^. Based 
on the magnitude of the effects in the DMT model and 
particularly in the SL model jBl) , we used a very simple 
double dipole approximation for the pion cloud 



Model qi, q 2 A^, Al iVi,T(0) G* M (0)/3, X 2 



G7T 
M_ 

3Gn 



A™ 



A: 



A 2 +Q 2 



(6.13) 



where G D = 1/(1 + Q 2 /0.71) 2 (with Q 2 in GeV 2 ) and A w 
and A 2 are parameters to be adjusted to the data. The 
parameter A^ can be interpreted as the fraction of pion 
cloud effects at Q 2 = and A 2 measures the falloff of the 
pion cloud. Note that we parametrize the ratio of G^ to 
3Gd following the tradition of scaling G* M with the dipole 
factor Gd- This form assume a falloff of G\ { ~ 1/Q 8 , to 
be compared with G* M ~ 3Gd 



1(B) 


0.169 




2.88 


0.548 




0.489 




0.792 


1.32 


11(B) 


0.181 




3.05 


0.547 




0.493 




0.790 


1.26 


I 


0.313 


0.474 


2.88 


1.026 




0.374 


1.172 


0.798 


2.64 


II 


0.290 


0.464 


2.95 


1.012 




0.393 


1.224 


0.794 


1.84 



TABLE II: The dimensionless Delta wave function parameters 
Qfi and OL2 and the normalization constant N\ are denned 
in Eq. ([27T6J. Model 1(B) and 11(B) include no pion cloud; 
Models I and II include a pion cloud with parameters A^- and 
Al (in GeV 2 ) defined in Eq. (|6TT3)l . The overlap integral 
between nucleon and Delta wave-function (which cannot be 
larger than 1) is X(0); G* M (0)/3 measures the quality of the 
fit for Q 2 = (where the experimental result ~ 1). 




i/Q 4 



FIG. 4: Fit to the G* M and Bare data using using model I 
and II. The Bare data is from [ill . l52l ]. G* M data from figure 
[3] The Bare result is now defined considering A,r = in the 
dressed models. 



It is important to realize that, due to the complexity 
of the pion production process, the decomposition of Eq. 
(|6.12p is strongly model dependent. Dynamical models 
can differ in coupling constants, off-mass-shell extrapo- 
lations, and off course the parametrization of bare com- 
ponent itself. Using Eq. (|6 . 1 3[) we can find several com- 
binations {Gm,Gm) with approximately the same sum 
G* M . We must find some way to constrain one of these 
two components. 

To do this, we use a procedure implemented for the 
first time in Ref. [4l[. Using the SL model, Julia-Diaz, 
Lee, Sato and Smith extract, independently at each Q 2 
point, a value of the bare form factor. This is possible be- 
cause the bare form factor, G^, is one of the parameters 
that enters the dynamical SL model, and it is therefore 
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possible to determine it, without any theoretical bias, by 
a best fit to the data. 

The data for the bare component of the form factor, 
determined in this way, are shown in Fig. 2J Note that a 
separation between the " bare" data and the experimental 
data can only be made in the region Q 2 < 2.5 GeV 2 . 
Above Q 2 ~ 3 GeV 2 the pion cloud contributions are 
less significant and the dynamical model produces only a 
correction to the bare component; in this region the full 
contribution comes mainly from the bare component. 

Models I and II, dressed by the pion cloud, are a result 
of a simultaneous fit of both components of the form 
factor Eq. (|6.12p , to the G* M experimental data [HI [H, 
l53l . [54J and of the bare component, to the bare "data" 
for Q 2 < 2 GeV 2 based on the SL extraction [U H| 
discussed above. The parameters for Models I and II arc 
compared to models 1(B) and 11(B) in Table HH (The x 2 
given there is for the fit to G* M only.) In Fig. [3] we present 
the results our predictions for the dressed models and the 
respective Bare version obtained setting A w = 0. In the 
same figure we can see that the data below 0.13 GeV 2 
(first three points) cannot be described by our model. 
This limitation is related to the behavior of the overlap 
integral of the nucleon and Delta wave functions. Note 
that the values of 1(0) are almost identical for the bare 
and dressed models, but the dressed models now have 
G* M (0)/3 « 1 due to the addition of a pion cloud term 
of about 46% at Q 2 = 0. Furthermore, a reasonable 
description of the 'bare data' is obtained for both models 
(x 2 = 4.2 for model I and \ 2 = 4.6 for model II) at least 
for Q 2 > 0.13 GeV 2 . 

It is worth to mentioning that the parametrization of 
the dressed models changes the results of the bare contri- 
bution relative to the Fig. [3J Although similar, the bare 
results presented in Fig. Q] are slightly larger than the 
results of the models 1(B) and 11(B) presented in Fig. O 
in particular for Q 2 < 2 GeV 2 . This increment is not 
obvious in the graphs, but as we can see in tabic HQ the 
parameters ct\ and «2 for the models 1(B) and 11(B) are 
significantly different from the models I and II. This fea- 
ture is the result of including the low momentum "bare 
data" (Q 2 < 2 GeV 2 ) that appears not be completely 
consistent with the high Q 2 data (Q 2 > 3 GeV 2 ) at least 
with the naive parametrization (|6 . 1 3[) . To sort out this 
situation more high Q 2 data of high quality are needed 
(the current data set includes only 6 data points with 
Q 2 > 3 GeV 2 , compared with 26 data points for G* M and 
21 data points for Gf, for Q 2 < 3 GeV 2 ). 



B. Comparing with other models for the Bare form 
factors 



As mentioned above, the dynamical models need a 
phenomenological parametrization of the "bare" vertex. 
There is some freedom in the choice of the "bare" form 
factor, but the constraints of the quark models are usu- 
ally taken into account. This parametrization can be 




FIG. 5: Comparing Bare form factor parametrization with 
the 'bare data' from Refs. |4l|. [52| with models I and II, SL 
model [II and DMT modelggP 



done for each transition multipole G* (a — A1,E,C). 
For the SL and DMT models the "bare" form factors can 
be written 



G^ = Gf / (0)(l + aQ 2 )exp(-6Q 2 )/, 



(6.14) 



where G£f(0) < 2.07 fixes the contribution of the quark 
core for Q 2 = 0, a, b are positive parameters, and / 

for the SL model and / = 



1 



1 



for the DMT 



(M+m) 2 

model. For SL Gf f (0) = 2; for DMT Gt(0) = 1.65. All 
other parameters can be found in Refs. [48l . |43 | . 

The structure of the Utrecht-Ohio model [5(| for the 
bare form factor is incompatible with a pion cloud which 
is not zero for Q 2 = 0, and for this reason a direct and 
simple estimation of the bare form factor based in Eqs. 
(|6.12[) and (|6.13[) is not possible. 

The bare form factors use by of SL and DMT [obtained 
from the analytical expression (|6.14|) ] and our numerical 
result for models I and II are compared in Fig. \5\ We 
see that the bare SL form factor overestimates the data 
for Q 2 > 4 GeV 2 and that DMT always underestimates 
the data (suggesting a significant pion cloud contribution 
even for Q 2 > 3 GeV 2 , since for 4 GeV 2 the effect is still 
about 10%). We need to point out that the particular 
parametrization of both SL and DMT, in particular the 
SL model, was done before the data Q 2 
came available 



12] 



> 4 GeV 2 be- 
The new data showed the limitation 
of the particular parametrization, and was one of the mo- 
tivations for work presented in Ref. [4l[ , where the Bare 
form factor are adjusted for each Q 2 point. Note that 
our results are very similar to SL for Q 2 > 3 GeV 2 . 
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C. Comparing Model I with Model II 

In this work we consider two different models for the 
nucleon wave function as presented in Ref. @. We con- 
clude that the model II is the best model for the 7 AT — > A 
transition because is the one that better describes the 
high Q 2 data (meaning the pion cloud contribution is 
better described by the dipole form). Model I gives a 
slightly worse description of the data, but both models 
are almost indiscernible in the region Q 2 < 3 GeV 2 . This 
result is very interesting because the models are funda- 
mentally different in the description of the nucleon form 
factors. In the model I the isospin symmetry is exactly 
imposed with /i + = /i_ (Eq. (14. 5p with c + — c_), lead- 
ing to the failure of the description of the neutron elec- 
trical form factor (see Ref. [5|) and a consequent high x 2 
penalization as shown in table Q] As the — > A tran- 
sition is independent of /i+ (only isovector form factors 
contribute) a reasonable description of the form factors 
can be obtained for both models. 



D. Discussion 

Since our model includes only S'-waves in the nucleon 
and A wave functions, the comparison of our results to 
all models and frameworks used in the description of the 
electromagnetic N — * A transition has to be done with 
care. As mentioned above, we can only predict non-zero 
contributions to the dominant form factor G* M (55% of 
it for Q 2 = 0), while the experiments reveal two more 
non- vanishing, albeit small, form factors (G* E and G* c ). 
The limitation of our model is visible in its failure at 
high Q , where G* E is comparable with — G* M , according 
to pQCD. That regime is however out of reach of the 
present state-of-the-art measurements. 

On the other hand, our results are hardly compara- 
ble with the low momentum Effective Field Theory and 
Chiral Perturbation Theory [H HE H3, S H • ThoSe 
models include pion degrees of freedom consistently at 
least at one pion loop level. But the range of the predic- 
tion is limited to Q 2 < 0.25 GeV 2 , due to the expansion 
in terms of the terms of the small variables (pion mo- 
mentum, pion mass, difference of A and nucleon mass) 
. Besides the range limitation, the bare contribution is 
adjusted to the data around Q 2 — and is not really 
predicted from quark structure. This leads to a bare con- 
tribution significantly different from quark models. Note 
that Effective Chiral Perturbation Theory relies on an en- 
ergy scale parameter, usually A ~ 1 GeV, decoupling the 
short range physics (bare) from the long r ang e physics 
where the pion cloud is included. In Ref. 15 7| the pion 
cloud gives a positive contribution; in Ref. [59( the con- 
tribution is negative (bare contribution Gf / (0) = 4.04). 

As for lattice QCD data, the comparison is not yet 
conclusive: the (quenched) lattice data of Alexandrou et 
al [icj overpredict G* M in the chiral limit for Q 2 > 0.1 
GeV 2 . In principle, the quenched lattice data should be 



comparable to the bare form factor results, and hence 
it might be expected to be larger that the experimen- 
tal data, but instead it underestimates the data. It is 
not known yet either this discrepancy is due to the lim- 
itations of the quenched data or to the extrapolation to 
the physical region. The available results from full QCD 
(unquenched) are not adequate for an extrapolation [6l| . 

The soundest comparison to be made, then, is to other 
valence quark models. This is why in the previous subsec- 
tion we compared the magnitude of the Bare form factor 
obtained by us to the results from dynamical models and 
constituent quark models based on S'-wave wave func- 
tions. We are left then with comparing our results to 
the predictions of Light-Cone Sum Rules of Braun et al 
[42I ]. This formalism divides the main contributions of 
the form factors into two components: the soft contribu- 
tion falling with l/Q 6 and the hard contribution due to 
pQCD with a l/Q 4 falloff. The soft contributions, domi- 
nant in the intermediate region, are explicitly evaluated, 
using the nucleon asymptotic amplitudes (valence quark 
distributions), and an adjustable momentum range pa- 
rameter (Borel parameter). Their results describe G* M in 
the 3 — 6 GeV 2 region, but fail in the region of low and 
high Q 2 (where they have an almost constant slope) . Also 
they underestimate G* M at low Q 2 like the constituent 
quark models (at Q 2 — 1 GeV 2 the prediction for G* M 
only takes account of 60% of the experimental value). 

The success of our model in the description of the bare 
form factors (high Q 2 region) is related to our choice 
for the form of ipA m Eq. (|2.16[) . In particular, the extra 
power in the ct2+XA factor, comparatively to the nucleon 
wave function, plays a key role. It is also interesting to 
note that our results would be very similar to the results 
of the Ref. [42j had we used a A scalar function of the 
nucleon- type (see Eq. (|2.9p 1. This generates G* M with an 
almost constant slope and G* m /(3Gd) — 0.6. As men- 
tioned, our parametrization for the A scalar form factor 
ipA is consistent with G* M ~ l/Q 4 for large Q 2 , i.e., the 
pQCD prediction. For the AA form factors (a A in the 
initial and final state) our prediction, based on Eq. (|2.16p . 
has a l/Q 6 falloff for the dominant form factors at high 
Q 2 , instead of the l/Q 4 advocated by pQCD. We can 
find a good compromise for both the low Q 2 description 
and the expected pQCD behavior, by considering 



(a[ + XA)(a' 2 + Xa) (ai + xa)(«2 + xa) 2 ' 

(6.15) 

where cJ i are new range parameters, and a, b coefficients 
which balance the two regimes: a the deep Q 2 asymptotic 
region and b the low-intermediate Q 2 region. This com- 
promise requires b » a pushing the pQCD dominance 
to very far away. 

VII. SUMMARY AND CONCLUSIONS 

We have developed a systematic formalism for the 
description of baryon wave functions, built upon con- 
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stituent quark (flavor, spin and isospin) and baryon (spin, 
isospin) effective properties. The form of the nucleon 
and A wave functions with S'-wave orbital angular mo- 
mentum components only is presented. The formalism is 
manifestly covariant. The wave functions are covariant 
and transform like Dirac and Rarita-Schwinger spinors. 
The matrix elements are covariant, with the same form 
in all frames. 

One of our models (model II) describes the nucleon 
form factor data |5( and the dominant contribution to the 
N — > A electromagnetic transition. The results for G* M 
show a reasonable agreement with the data, and explain 
its measured falloff. They are also consistent with the 
long range behavior predicted by pQCD (G* M ~ l/Q 4 ) 
[32| . Our results for G* M are consistent with the results 
of quark models where only S'-states are considered. 

In agreement with previous works (see Refs. [TU, Hi], 
HUE!])) we conclude that a successful description of G* M 
near Q 2 = requires an addition of a pion cloud term 
not included in the class of valence quarks to explain 
the strength at Q 2 = of the magnetic form factor G* M 
of the TV — > A electromagnetic transition. Our predic- 
tions for the pion cloud underestimate the predictions 
based on the Sato and Lee model (46% versus 33%) for 
Q 2 = 0. This gap can in principle decrease once D-states 
are included into the A wave function. The magnitude 
of the pion cloud for Q 2 — is similar to the estima- 
tions of the DMT model but we predict a faster falloff. 
[Note that our results are not directly compared with the 
DMT model because our 'bare' contribution is fixed by 
the 'bare data' extraction of the SL model.] Except for 
the region Q 2 ~ 0, our model is consistent with the 'Bare 
data' extraction based in Sato and Lee model [4l|. For 
the region Q 2 > 3 GeV 2 our model and the original SL 
model 0| overestimate the data slightly. In this region 
either the pion cloud parametrization is not adequate or 
the data is insufficient to constraint adequately the pa- 
rameters of the pion cloud. More higher Q 2 data and 
more accurate data for both G* M and Gfj (SL model) is 
therefore necessary for to establish the effect of the pion 
cloud. 

Next, we plan to generalize the structure of the wave 
functions to include higher orbital angular momentum 
states in the quark-diquark system, without loss of the 
covariance requirement. The inclusion of D states in the 
nucleon and A is in progress. 
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APPENDIX A: NUCLEON WAVE FUNCTION 

In the non-relativistic limit, the wave function of a spin 
1/2 system composed of a quark (spin 1/2) and a quark 
pair (diquark), as in the nucleon, can be decomposed into 
two components: the scalar part (spin-0 diquark) and the 
axial- vector part (spin-1 diquark). 

The spin-0 part is just 

4>l=Xs, (Al) 

where Xs is the usual Pauli spinor (2x1 state) . The spin 
of the system is then given by the spin of the quark. The 
relativistic generalization is u(Po, s), where Po = {m, 0} 
is the four-momentum of a nucleon at rest. 

The spin-1 component of the wave function 4>\ de- 
scribes the quark-diquark spin 1/2 system (nucleon) in 
the initial state and a diquark polarization vector in the 
final state. The quark-diquark spin state, represented by 
Vn s , is a direct product of a spin-1 diquark state with a 
spin 1/2 quark state. In this case the three-component 
vector s\ (with i = x,y,z and A = 0, ± the diquark po- 
larization index) describes the diquark polarization [these 
vectors are the three-component parts of the polarization 
vectors defined in Eqs. (|3.1[) ] . Then 

{V Ns y = Y,^y\WxXs>, (A2) 

As' 

where s — ±1/2 is the spin projection of the nucleon, 
(si mi; s\ m2\j nij) is the Clebsch-Gordan coefficient that 
couples spins si and S2 to total spin j, and \s' the quark 
spinor. Explicitly 




Equation (|A2j) can be written 

(V Ns y = -\°iXs, (A3) 

where Xs now represents the nucleon spinor. The natural 
relativistic generalization is 

(V Ns y - U a (P Q ,s) = i= 75 (V - ^) u(P ,s), 

(A4) 

where u(Po,s) is the Dirac spinor of the nucleon and 
a = {0, i} with U° = 0. 

Note that Eq. (|A4|) describes only the initial state 
of the vertex represented in Fig. [T] (the 3-quark bound 
state). To obtain the amplitude of the full process we 
need to contract (V,) 1 with the diquark final state e"* 
(the quark is off-shell with its final spin state unspeci- 
fied). As result, we have the amplitude 

cfi =-eT Po U a {P ,s). (A5) 

Equations (|A4[) and (|A5[) can both be generalized for a 
moving nucleon by means of a boost in the z-direction. 
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APPENDIX B: SPIN STRUCTURE FOR THE A 
S-STATE 

In close analogy with the nucleon, the A wave function, 
in its rest frame rest frame, can be written as a direct 
product of a spin-1 diquark and a spin-1/2 quark 



(w = E< 1A ^'i§^ 



\Xs> 



(Bl) 



As' 



where s = ±3/2 or ±1/2 is the spin projection of the A. 
Once again e\ is the diquark polarization in the A rest 
frame, and Xs' a quark Pauli spinor. 

We can also express Va s in terms of a basis of spin 3/2 
states: 














1 





Ui i = 

1" 2 





W 




W 






(°\ 










1 


LU 3 = 
2 





W 




w 



(B2) 



In this case the connection between the u> and Va s can 
be written 



(v As y = 'Pus, 



(B3) 



where T l is an 2 x 4 matrix that transforms the spin 3/2 
state of the A into the spin 1/2 state (of a quark). The 
elements of T l can be evaluated using the coefficients in 
Eq. ifHTjl . The result is 



(B4) 



j_/\/3 -1 

'y/6 \ 1 

J_fV3 1 
VQ\0 1 V3y 

2/0 1 0' 
3 I 10, 



As in Eq. (|B1|) . (Vas) 1 is a 2 x 1 spinor with a spin 1/2 
structure. 

To convert (|B1|) to relativistic form, add a negative 
energy 4x1 lower component that vanishes in the A rest 
frame: 



r z 




(VAsT 







w\P ,s). 



(B5) 



Here w l (Po,s) is the Rarita-Schwinger vector-spinor for 
3- momentum P = 0. It satisfies the constraint Eqs. 
(|2.12|) . These constraints insure that w° vanishes in the 
rest frame, as implied by Eq. (|B5[) . 



To generalize the states to an arbitrary frame with 
P / 0, wc boost them using a Lorentz transformation A, 
giving 



■/(P, S )=S(A)AP a w a (P Q ,s) 



(B6) 



Using the state w@ in a arbitrary frame, the A-quark- 
diquark vertex is constructed in the same way as the nu- 
cleon vertex. Considering the final state diquark polar- 
ization vector, e^, following the nucleon state convention 
of Eq. C£2J| 



(B7) 



gives the spin wave function introduced of Eq. (|2.11[) . 



APPENDIX C: RELATIVISTIC SPIN 
PROJECTION OPERATORS 

In this paper we work with operators O a ^ that satisfy 
the constraint equations 

P a O aP = = O aP P p (CI) 

In the particle rest system, such operators "live" in the 
3x3 subspace corresponding to nonrelativistic 3 dimen- 
sional space, and it is easy to relate these operators to 
their nonrelativistic analogues. 

As an example, consider the projection operators that 
operate on the direct product of spin-1 and spin-1/2 
spaces. The total angular momentum operator is the 
sum 



r = w i + s i , 



(C2) 



where W l are spin-1 operators (with multiplication by 
the unit operator on the spin- 1/2 space implied) and S l 
are the spin-1/2 operators (with multiplication by the 
unit operator on the spin-1 space implied). The projec- 
tion operators are constructed from the operator 



2W- S 



W 2 



s 2 



1 J 

-2 J 



(C3) 



Hence the projection operators Vj are 
P 1/a = |(l-2W-S) 



V. 



3/2 



Using (W jk ) 1 = -ie ljk and S l 



2W-S) . 

= <t 1 /2 we get 



(C4) 



(Pz/2)jk = S jk - | (Jj<Jk , (C5) 
leading immediately to the relativistic generalizations 

I'D V* l ( FoPoV ( P P \ 

(V 3 /2) a - 9 a p ~ - (Vi/2) a p ■ (C6) 
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Boosting these from Pq to P gives the operators intro- 
duces in Eq. (|2.23l) above. 

References [26J, |27| define two spin 1/2 projection op- 
erators. In addition to V-yi^ (which they call V\\) they 
introduce the operator 



pap 







M 2 



(C7) 



is a spin 1/2 projector on the (0,4) space. Since we 
work in the space of operators satisfying the constraint 
(|Cip . this operator is excluded from our basis. In our 
formalism, it is part of the operator g. 

We conclude this appendix by proving the useful rela- 
tion (|4.16[) . Start by using the properties 



PtD„ 



= D^P V _ 



(C8) 



which follow directly from the definition of as a sum 
over polarizations. Then we prove two intermediate re- 
sults. First, using the notation V^? 2 (P+) to denote the 
projection operator with momentum P + for a state with 
mass M [and similarly for P?Y 2 (P-)], we see that 

V?; 2 (P + )D a pV^ 2 = IT(P + ) [i a D a ^]r{P-) . (C9) 
If b = P + ■ P_ , the quantity in square brackets is 

b ^ 



[■■■] 



-4 + 
-4 + 



b 

b ~ 

M 



m 



-4 - 2ab + 



ab 2 
Mm 



(CIO) 



where we used the fact that — > — to when acting 
to the right, and f + — > — M when acting to the left, 
because both anticommute with the 7 standing to the 
right and left, and then can be eliminated using the Dirac 
equation satisfied by the incoming and outgoing states. 
Then, using f_f+ = 26 - -> 2b - Mm, and the 

value of a from Eq. (|4.10[) , we get 

[•■■]=-3 (Cll) 
and Eq. (|C9|) reduces to 

V?; 2 (P + )D a pP^ 2 = -lr(P + )T(P-) ■ (C12) 
Similarly, using the same procedure we can show that 

gTD a pPlJ 2 = -lY(P+W{P-) ■ (Ci3) 



Combining these results gives the result we seek 

-V?; 2 (P + )D al3 V^ 2 = 0. (C14) 



1 3/2\ r +l lya P l 1/2 — y 



APPENDIX D: AN EXAMPLE OF THE 
IMPORTANCE OF THE COLLINEARITY 
CONDITION 

To give some insight into the importance of the 
collinear frame in the definitions of fixed-axis polariza- 
tion states, consider a simple example where the initial 
and final baryon are identical (both have mass m) and 
the four-momenta are not collinear 



P± = (E',psm6,0,±pcost 



(Dl) 



with E' = \J m 2 + p 2 . We can transform these momenta 
to a collinear frame by boosting in the x direction using 
the transformation 



D , 



cosh?7 sinhf? 

sinh?7 cosh rj 

10 

1 



(D2) 



Collinearity is achieved if 

sinh r\ E + cosh rj p sin 9 = 0, 
and the resulting collinear four momenta are 
P ± =B X P' ± = (E,0,0,±pcos9) 



(D3) 



(D4) 



with E = ^Jm 2 + p 2 cos 2 9. 



In this example, consider the longitudinal polarization 
vectors only. In the collinear frame they are 



1 



-'0F± 



(±pcos6»,0,0,E) . 



(D5) 



Their scalar product is 

e op = -p 2 cos 2 9 - E 2 



m 2 s n 



P+P-. (D6) 



In the original, non-collinear frame, the longitudinal po- 
larizations are 



£ op' - B x £ OP ± 



(± cosh 77 p cos 9, =p sinh r\ p cos 9,0, E) . (D7) 



Hence, 



2 

me , ■ e . 



= -(p 2 cos 2 9 + E 2 ) 

= -(2p 2 cos 2 9 + to 2 ) = -Pi • P'_ (D8) 



showing that the scalar product of the two longitudinal 



polarization vectors e 



is invariant and uniquely 



defined. 

Now, suppose we were to define the longitudinal vec- 
tors in the original frame. One way to do this would be 
to observe that the two momenta (|D1[) can be obtained 
from the vectors 



P± = {E',0,0,±p) 



(D9) 
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by rotation about the y axis by an angle ±9 



l 

cosf 






±sin( 

1 



(D1Q) 



=F sin 9 cos 9 
The longitudinal vectors corresponding to (|D9[) are 



1 



(Dll) 



m 

so the new polarization vectors would be 
e ± = R±e e± 

= — (±p,± sin 0£J',O, cos <?£'). (D12) 
m 

These vectors are completely different than the correct 
ones given in Eq. (|D7[) . and in particular 

to 2 £'+■?'_ - -(p 2 ~sin 2 9E' 2 +cos 2 9E' 2 ) 
= -(2p 2 cos 2 9 + m 2 cos 20) 
^ -P' + -P'_ (D13) 

In conclusion: the correct way to treat fixed-axis po- 
larization vectors is to transform to a collinear frame (if 
necessary), define the fixed-axis vectors there, and then 
do the inverse transformation back to the original frame 
(if desired). 

APPENDIX E: IN THE COLLINEAR FRAME 



Consider the tensor 



e XP + £ XP- 



(El) 



where both e\p + are polarization vectors initially ori- 
ented along the z axis (in a collinear frame) and satisfy 
the constraints P+ • ep + = P- ■ ep_ = with 

P± = (E±,0,0, P ±) (E2) 

and P 2 = M 2 , P 2 = m 2 . is a sum of direct products 
of four-vectors, and therefore is a covariant tensor. The 
explicit form of the £\p ± is 

±=(0,-A,-i,0) if A = ± 





l 




1 


£op_ 


= — 




m 




1 


£op + 


~ M 



(p+, 0,0,^4 



(E3) 



Using these explicit forms and interpreting e p as a 
column vector and e P *_ as a row vector, we get the fol- 
lowing matrix form for D 



Mm 



p+P- P+E_ 
Mm 
Mm 

E+p- E+E- 



(E4) 



The covariant form for this tensor can be found by 
exploiting the fact that P+^D^ = and D» V P- V = 0. 
Hence the most general form of D^ v is 



ax -<r + 



p^Pi 



where b = P + ■ P_ and a\ — 1 (to give the correct D xx 
and D yv components). The coefficient can be found 
from the trace 

D\ = -2- ?±-^- 

M Mm 



which gives 



Mm 



« 2 



b(Mm + b)' 



It is easy to verify that the two forms (|E4j) and (|E5|) are 
identical. 



APPENDIX F: THE CURRENT J" FOR THE 
7 N^ A TRANSITION 



Equations (|5.18|) and (|5.19[) can be written 

'AN\ 



1 



--^fvWp(P+,X+) 



where the form factor 



7 ^ q 7q7 5 J U (P_,A_)X, 

(Fl) 

was defined in Eq. (|5.2[) and 



I 



i> A (P+,k)iP N (P-,k). 



(F2) 



The operator in the square brackets in Eq. (|F1[) is re- 
duced using the form of D° a given in Eq. (|4.9p . remem- 
bering that the properties of the Rarita-Schwinger wave 
function imply that terms proportional to P+ and 7^ 
(when operating to the left) are zero, and that the Dirac 
equation may be used to replace f _ — * m when oper- 
ating to the right, and f + — ► M when operating to the 
left. We get 



V- V(^+-m) 

2A(M^ - P£)/ - 2g»P 



where 



A = 



Mm + b (M + m) 2 + Q 2 ' 



(F3) 
(F4) 



Noting that P^ = P>* + , the operator (jF3| can be 
written in terms of the invariants of (|5.11[) 

$i«V + Wi* + WV - <?4<?H 7 5 (F5) 
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where 



.9! = 2MA 

92 = -2A 

93 = -A 

94 = 2. 



Note that 



,g 4 = (M + m)g\ + 



M 2 - m 2 



-92 - Q 2 93 



(F6) 



(F7) 



as required by current conservation, Eq. (|5.13p . The 
Jones and Scadron form factors are 



Gi — fv9i 



y/3' 



(F8) 



We conclude that the physical form factors are, within 
the 5*-wave model, 



8m 



G * MiQ2) -3V3(M + m) 
G* E (Q 2 ) = G* c (Q 2 ) = 0- 



(F9) 
(F10) 



APPENDIX G: ASYMPTOTIC Q 2 DEPENDENCE 
OF THE INVARIANT BODY INTEGRALS 



In this appendix we discuss the asymptotic dependence 
of the "body" integrals 



B H (Q 2 



tf> H (P + ,k)il> N (P-,k), 



(Gl) 



where H = N or A. [5[. The high Q 2 dependence of 
Bh(Q 2 ) determines the asymptotic behavior of the nu- 
cleon and N — > A form factors. To simplify the dis- 
cussion we consider the easiest case, when the parame- 
ters of wave function are Pi = @2 = P for the nucleon, 
a\ = oii = a for the A, and we will sometimes use the 
notation Ph, where Pn = p and Pa = ct. 

The integral B is covariant and may be evaluated in 
any frame. It is convenient to evaluate it in the "anti- 
lab" frame, where the final hadron is at rest. In this 
case the momenta are Pjj = {mn, 0,0,0) and P qo = 
(Eq, 0, 0, — go)j with Eq = \J m 2 + q 2 , the nucleon energy 
in the initial state. The photon four-momentum is then 
q = (m H - E o ,0,0,q Q ) with 



V 2 o = 



Q 2 



Q 4 

Am 2 H 



(G2) 



as Q 2 —> oo. With these momenta we can write the body 
integral as 



B H (Q 2 ) = 



N 
(2*0 



k 2 dk 
2m„ E> 



^ H (P H ,k)I(Q 2 ), (G3) 



where 

HQ 2 ) = 



dz 



" mm. m. m.„ 



(G4) 



In this frame only the initial state depends of the angular 
coordinate z = cos#. Introducing the parameter 

(J3- 2)mm s + 2E E S 



gives 



HQ 2 



Aq 2 k 2 
™2 



2q G k 



dz 



(z + V ) 2 2q 2 k 2 {l 1 2 -l) 



2% 2 



(G5) 



as Q 2 — > oo. Motivated by the nonrelativistic definition 
of the wave function at the origin, we define the following 
covariant integral 



<M0) = 



k 2 dk 
2m , E, 



ii H (P H ,k). 



(G6) 



Our results can now be expressed in terms of the behavior 
of this integral. 

Case I: If the integral (|G6[) exists, then the diquark 
momentum k is localized, and the limit (|G5p can be taken 
under the integral, leading to the result 



lim B H {Q 2 ) 



Nq 
(27T) 2 
Nq 

(2tt) 2 



^h(0)I(Q 2 ) 
^(0)— ^— . (G7) 



We obtain the interesting (and well known) result that, if 
cases where the value of the relativistic wave function of 
one of the hadrons is finite at the origin, the asymptotic 
from factor is determined by the high momentum behav- 
ior of the other wave function. For the models used in 
this paper this shows that the N — > A body form factors 
go like Q~ 4 are large Q 2 . 

Case II: If -0(0) does not exist, the limit (|G5[) cannot 
be taken and the analysis of the large Q 2 behavior de- 
pends on the behavior of the full integral. In this case 
we return to (|G3[) and write (H = N now) 



(G8) 



The integral X can be evaluated in the diquark energy 
scaled by the diquark mass x = E s /m s considering k = 
m s \Jx 2 — 1. As result we have for H = N 

g( x ) 



AQ Z ) = 



! D(x + uj) [ 



;dx, 



(G9) 
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where g(x) = V x 2 — 1; lu = 2) and the denominator 
D is 



D = 



in 



rn- 



X + ZLO X H r. 

m m z 



(G10) 



The last approximation holds at large go, where all con- 
stant terms and terms proportional to x/qo can be ne- 
glected, because they are always smaller than terms pro- 
portional to x 2 , (qo/m) 2 , or xqo/m. Also for large qo we 
can re- write (IG9I) as 




Note that each term inside the brackets diverges but the 
result is convergent. The above integrals can be per- 
formed analytically following the usual techniques. Con- 
sidering x = 1/ cosu the integrand function becomes an 
algebraic function of cos u and sin u that we integrate an- 
alytically using the Mathematica program. Considering 
only the leading and next leading terms in qo/m in the 



the general expressions, one is left with 



,2 r 



?o 2 L 



( — ) Ru-) 
1 m ' 



where 



n(w) = i + 



(Gil) 



(G12) 



2 tan" 1 



1 -LU 



tan 



The analytical continuation of the (|G12[) for the case 
u> > 1 (or (3 > 2) is obtained considering the relation be- 
tween logarithms and arc-tangent log j^ff = 2i tan _1 (.T) 



and the replacement yl — lu 2 
In conclusion, we can write 



B N {Q 2 



N 2 



lo 6 3 



(An) 2 Q 4 b m 



(G13) 



This term is independent of (3. The nonleading terms in 
Q~ 4 carry the dependence in the parameter j3 and set 
the scale of the logarithm behavior. 

This logarithmic dependence of the nucleon form fac- 
tors was missed in Ref. [29], but this oversight does not 
affect any of the conclusions of that paper. 
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